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We investigate the dependence of elliptic flows vi on transverse momentum P T for charged hadrons produced 
in nucleus-nucleus collisions at high energy by using a multi-source ideal gas model which includes the 
interaction contribution of the emission sources. Our calculated results are approximately in agreement with 
the experimental data over a wider P T range from the STARand ALICE Collaborations. It is found that the 
expansion factor increases linearly with the impact parameter from most central (0-5%) to mid-peripheral 
(35-40%) collisions. 
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I. INTRODUCTION 

The first elliptic flow data from the Large Hadron Collider(LHC) obtained by the ALICE Collaboration fll are 
similar as observed at lower energies at the Relativistic Heavy Ion Collider (RHIC) flS]. The elliptic flow of 
charged hadrons as a function of transverse momentum V2(Pt) increases with Pj and saturates at higher Pj. The 
similarity between the elliptic flows at LHC and RHIC is consistent with the predictions of the viscous hydro- 
dynamic model |0]. The elliptic flow is generated in a collective expansion of the dense matter created in the 
collisions. The comparison to model calculations can provide valuable information of the collisional evolution. 
Recent results of relativistic viscous hydrodynamic calculations for Pb-Pb collisions at y/s^N - 2.76 shows sys- 
tematic deviations from the ALICE data in the small momenta region Pj < 800 MeV 15[]. The deviations is 
interpreted to be a non-negligible contribution of non-thermalized particles from jet fragmentation. The origin of 
this deviation is still under debate. Elliptic flow results from the interactions among the produced particles and can 
be used to probe local thermodynamic equilibrium. By fitting the experimental data of elliptic flows, it has been 
found in Ref. |Q] that the average Quark Gluon Plasma (QGP) specific shear viscosity rj/s slightly increases from 
RHIC to LHC energies in the newly developed hybrid model VISHNU, which connects viscous hydrodynamics 
with a hadron cascade model. At the LHC energies, it remains an open question whether the QGP fluid is more 
viscous or more perfect. 

Elliptic flow is a measure of the azimuthal anisotropy of particle momentum distributions in the plane 
perpendicular to the beam direction. It results from the initial spatial anisotropy in non-central collisions 
and is thus sensitive to the properties of the hot dense matter formed during the initial stage of heavy ion 
collisions. In a thermodynamical picture, the asymmetric distribution of the initial energy density causes a 
larger pressure gradient in the shortest direction of the ellipsoidal medium, which can lead to an anisotropic 
azimuthal emission of particles. As a signature of collective flow in relativistic nuclear collisions, the elliptic 
flow was experimentally observed at the AGS, SPS, RHIC, and LHC,respectively Jz-Ll]. In recent years, it has 



been believed that the system created at the RHIC 111 10 is a strongly coupled quark-gluon plasma (sQGP) and 
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behaves like a nearly "perfect" fluid by using the hydrodynamic simulation to the elliptic flow at the RHIC 

II. MULTI-SOURCE IDEAL GAS MODEL 



The multi-source ideal gas model is developed by us Ill4l4l6ll from the Maxwell's ideal gas model. In this 
section, we briefly outline the calculation in the model. Many primary nucleon-nucleon collisions happen at the 
initial stage of nucleus-nucleus collisions. Each primary nucleon-nucleon collision can be regarded as an emission 
source (a compound hadron fireball) at intermediate energy or a few sources (wounded partons and woundless 
partons) at high energy. The participant nucleons that take part in primary collisions have probabilities to take 
part in cascade collisions with latter nucleons. Meanwhile, the particles produced in primary or cascade nucleon- 
nucleon collisions have probabilities to take part in secondary collisions with latter nucleons and other particles. 
Each cascade (or secondary) collision is also regarded as an emission source or a few sources. Many emission 
sources of final-state particles are expected to be formed in the collision process, and each source is considered 
to be a Maxwells ideal gas in the rest source frame. There are interactions among the emission sources due to 
the mechanics and electromagnetism effects. The interactions will cause the source to depart from the isotropic 
emission. To describe the anisotropic source, the sources deformation and movement in the model have been 
introduced. 

Let the beam direction be the oz axis and the reaction plane be the xoz plane, which is spanned by the vector 
of the impact parameter and the beam direction. In the rest source frame, the three momentum components P' x , 
P' and P' of the considered particle have Gaussian distributions with the same standard deviation cr. The az- 
imuthal anisotropy of particles produced in nuclear collisions can be described by the two momentum, P x and P y 
components in the transverse plane. The P x and P y distributions are given by 



fp-JP' x , y ) = -7=^exp 
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According to Gaussian momentum distribution, the width of the distribution cr = yTM, where T and M are the 
source temperature and particle mass respectively. 

Due to the interactions with other emission sources, the considered source will have deformations and move- 
ments along the ox and oy axes, we have 



P x = A X P \ + B x , (2a) 
P y =A y P' y +B y , (2b) 

where A A v and B xy represent the coefficients of source deformations along the ox(or oy) axis and movements 
along the corresponding axis respectively. The shifted deformable source is simply parameterized to have the 
linear relation between P xx and P' x v , which reflects the mean result of source interaction. For a Maxwells ideal 
gas, A XJ = 1 and B XY = 0, the source has no expansion and no contribution to the collective flow. For A v v > 1, the 
source expands, leading to an azimuthal anisotropy of the particle momentum distribution in the transverse plane. 
The amplitude of the azimuthal anisotropy can be measured by the elliptic flow, which is a sensitive probe of the 
final state interactions. For B xy , a positive value means source movements along the direction of the ox or oy axis, 
and a negative value indicates source movements along the opposite direction of the corresponding axis. 

Generally speaking, different A xy or B xy can be obtained for events with different centralities or impact param- 
eters. When two ultrarelativistic nuclei collide at non-zero impact parameter, their overlap area in the transverse 
plane has a short axis(ox), which is parallel to the impact parameter, and a long axis which is perpendicular to the 
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impact parameter. Elliptic flow in the central region of collisions is driven by the almond shape of the overlapping 
region I17H19I1 . This almond shape of the initial profile is converted by the pressure gradient into a momentum 



asymmetry. More particles are emitted along the short axis and the magnitude of this effect is characterized by A x . 
The momentum anisotropy is largest in the early evolution of the collision. As the system expanding, the system 
becomes more spherical, and the driving force quenches itself gradually. With the increasing of P x or P y of the 
particles produced in the collisions, the momentum anisotropy decreases. So, for the concerned centrality, A x is 
not a constant, but a function of the momentum. We set empirically the function 



A x (P x ,P y ) = 1 + £exp(- yl( Ul P x ) 2 + (m 2 ^) 2 /o-) 



(3) 



where the expansion factor k indicates the magnitude of the expansion, u\ and m 2 denote the contributions of P x 
and P, respectively. When k > 0, we have A x > 1, which renders the expansion of sources along the axis. When 
P x (or P y ) increases to sufficiently large values, A x closes to 1, which renders the momentum anisotropy vanishes 
and the elliptic flow v 2 = 0. The other parameters taken for calculations are A v = 1 and B x = B y — 0, which is the 
default value. 

The probability density function of P xy is obtained as 

</',, - />\,r 
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The corresponding transverse momentum and azimuth read as 
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respectively. With Eqs. (??), the joint probability density function of Pj and is given by 
fp T ,<p(PT,4>) = Pt/pt^Pt cos (p,P T sm(p) 
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Considering the measurement in experiments, the elliptic flow in the model is defined as 

v 2 (Pr) = <cos(20)> 



if d( f> fpr-LAPr ^s(2cf>)f PT ^P T , 4>)dP T 
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where APt is a given Pj bin. In the Monte Carlo calculation, we have 



P ' Xty = cry/-2 lnr u cos(27!r 2 , 4 ) (9) 

where n, r%, r^, and r\ denote random numbers in [0,1]. Using Eq. (2), (5), and (9), the transverse momentum and 
azimuthal angle can be represented as 

P T = yj[A x cr V-21nn cos(2^r 2 ) + B x j~ + [a v ct V^^lno cos(27rr 4 ) + B y j~ (10) 

and 



AyO" V-2 In cos(2nr4) + B y 

4> = arctan — -^ = ^ = , 

A K cr V-2 In rj cos(27rr2) + B x 



(ID 



4 



respectively. The elliptic flow reads 

\Pl-P 



III. COMPARISONS WITH EXPERIMENTAL DATA 

The dependences of elliptic flows on transverse momentum Pj for charged particles produced in Au-Au 
collisions at ysjvjv - 200 GeV are presented in Fig. 1 The symbols are the experimental data of the STAR Collabo- 
ration J3l in eleven collision centralities. Our results are shown as the curves. The parameters used for calculations 
and the corresponding^ 2 per degree of freedom (j^ 2 /dof) are given in Table I. In the calculation, the values of cr 
and k are determined by fitting the model to the data, the values of u\ and ui are taken to be 0.944 and 0.59-0.61 
respectively, which is independent on centrality. It is obvious that v'2 increases with Pj in the low Pt region, as 
predicted by ideal hydrodynamic calculations. The observed V2 saturates or decreases in the region of Pj > 2 



GeV/c. However, ideal hydrodynamic model calculations show that V2 increases with Pt in this region [20]. One 
can see that the calculated results are approximately in agreement with the experimental data on the dependence of 
V2 on Pt in the whole observed Pt region for all concerned centralities. The expansion factor k increases and the 
momentum distribution width cr decreases with the centrality percent. According to cr — ~\/TM, we can say that 
the source temperature T increases with the centrality. 

In Fig. 2, we show elliptic flows V2 as a function of the transverse momentum Pt for charged hadrons produced 
in Pb-Pb collisions at y/s^N = 2.76 TeV U. The symbols are the experimental results obtained with 4-particle cu- 
mulant methods and denoted as V2{4). The elliptic flow of charged hadrons measured by the ALICE Collaboration 
as a function of transverse momentum vi{Pt) is nearly identical to that measured by the STAR Collaboration at 
RHIC up to Pt — 3 GeV/c, independent of collision energy and centrality. The calculated results are approximately 
in agreement with the experimental data on the dependence of V2 on Pt for ft* in Pb-Pb collisions at yfsm = 2.76 
TeV. In the calculation, we take u\ = 0.988 and U2 = 0.625 - 0.640. The values of the expansion factor k increase 
with the centrality percent, and are systematically larger than those in Au-Au collisions at -\fstw = 200 GeV. The 
values cr and T also decreases with centrality percent, as that in Fig.l. 

The square values of the expansion factor k 2 used for Fig. 1 and Fig. 2 are displayed in Figs. 3, by different symbols 
as marked. In Au-Au and Pb-Pb collisions, Vz(Pt) tends to a saturation for the more peripheral collisions(> 40%). 
We find that, except for the saturation region, the values of k 2 exhibit a linear dependence on centrality percentc. 
For Au-Au and Pb-Pb collisions, we have k 2 = (3.161 ± 0.002)c and k 2 = (3.492 + 0.012)c, respectively. By the 
geometric relation of centrality percent c to the impact parameter b, we have c oc b 2 , which holds to a very high 
accuracy for all but most peripheral collisions. Thus, the expansion factor k linearly increases with the increasing 
of the impact parameter b from most central (0-5%) to mid-peripheral (35-40%) collisions. 

IV. CONCLUSIONS 

In the above discussions, we have investigated the results of azimuthal anisotropy by using the elliptic flow 
V2 vs the transverse momentum Pt for several centralities in y/Im = 200 GeV Au-Au and y/sm = 2.76 TeV 
Pb-Pb collisions. The calculated results are obtained in the framework of the multi-source idealgas model by 
comparing the experimental data. The calculated results are approximately in agreement with the experiment data 
of the STAR and ALICE Collaboration. In our model, the local sources of final observed particles are formed in 
heavy ion collisions. Their interaction are related to the hot dense matter in the sources, and also results in the 



azimuthally anisotropic expansion in the momentum space [21]. The almond shaped interaction volume produced 
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in a non-cental collision is converted by the pressure gradient into a momentum asymmetry. The parameter A x in 
the model is used to reflect the expansion of momenta, as shown in Eq.(3). An isotropic emission corresponds to 
k = 0. 

From the above comparisons, we can see that the values of the expansion factor k increases linearly with the 
impact parameter b increasing over the range of collision centrality (from 0-5% to 35-40%) in -^sJw = 200 GeV 
Au-Au and yfsm = 2.76 TeV Pb-Pb collisions. It supports the idea that the increase from central to peripheral 
collisions reflects the expected increase due to the change in initial eccentricity of the initial fireball from central to 
peripheral events at each centrality 1122.1 12311 . The difference between the two collisions is that the expansion factor 
of the latter one is larger than that of the former one. From RHIC to LHC energies, it is found that the average 
QGP specific shear viscosity tj/s slightly increases Recent hydrodynamic calculations in the framework of 

relativistic dissipative hydrodynamics have been interpreted as the evidence that the elliptic flow at RHIC energies 
is insensitive to the QGP viscosity and only becomes sensitive to it at LHC energies 12511 . 

The azimuthal anisotropy resulting from the final-state particles is one of the most informative quantities in better 
understanding the nature and properties of the matter in high energy nuclear collisions. The ideal hydrodynamic 
model calculations reproduce the mass ordering of V2 in the relatively low Pj region, but overshoot the values of 
V2 for all centrality bins 126J. To understand the viscous nature of QGP, dissipative hydrodynamics has recently 
been applied to explain the experimental data of the collective flow parameter v% by including the effect of shear 
and bulk viscosity I27M29I1 . Analysis of the elliptic flow in our simple model, where the system expansion can 
be quantified in the momentum space, shows that the expansion factor k is characterized by the impact parameter 
b, which is related to participating nucleons using a realistic description of the nuclear geometry in a Glauber 
calculation J30I1 . 

In summary, in the framework of the multi-source idealgas model, we have investigated the transverse mo- 
mentum dependence of elliptic flow for charged particles in Au-Au collisions at the RHIC energy and in Pb-Pb 
collisions at the LHC energy. Our calculated results obtained in the multi-source idealgas model approximately 
agree with the experimental data. The expansion factor k used in the calculation exhibits a linear dependence on 
the impact parameter b from most central (0-5%) to mid-peripheral (35-40%) collisions. Elliptic flow has been 
proved to be very valuable for understanding relativistic nuclear collisions. As we know that the hydrodynamic de- 
scription has gained important result over the past few years. The present work shows that the simpler multi-source 
idealgas model also can do the similar job. Particularly, in the descriptions of (pseudo)rapidity and multiplicity 
distributions for produced particles, this model is successful. In the description of the vi{Pt) at higher energies, 
the work is a successful attempt. 
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TABLE I: Values of parameters k and <x in our calculations. 



Figure 



Centrality cr(GeV/c) jfc(GeV/c) rVdof 



Fig.l 



Fia.2 



50%-60% 
45%-50% 
40%-45% 
35%-40% 
30%-35% 
25%-30% 
20%-25% 
15%-20% 
10%-15% 
5%-10% 
0%- 5% 
40%-50% 
30%-40% 
20%-30% 
10%-20% 



1.119 ±0.006 
1.127 ±0.004 
1.141 ±0.011 
1.164 ±0.016 
1.204 ±0.024 
1.262 ±0.025 
1.336 ±0.028 
1.365 ±0.020 
1.402 ±0.019 
1.452 ± 0.016 
1.551 ±0.016 
1.192 ±0.020 
1.248 ±0.022 
1.344 ±0.018 
1.486 ±0.015 



1.076 ±0.001 
1.069 ±0.002 
1.059 ±0.005 
1.046 ± 0.008 
0.991 ±0.010 
0.936 ±0.012 
0.875 ± 0.024 
0.719 ±0.052 
0.586 ± 0.045 
0.426 ± 0.042 
0.264 ± 0.036 
1.116 ±0.002 
1.110 ±0.004 
0.942 ± 0.056 
0.725 ± 0.044 



0.707 
0.716 
0.643 
0.622 
0.495 
0.413 
0.244 
0.192 
0.165 
0.146 
0.104 
0.808 
0.244 
0.124 
0.116 
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FIG. 1 : Pj dependence of V2 for charge hadrons in Au-Au collisions at ^js NN = 200 GeV. Experimental data taken from the 
STAR Collaboration U are shown with the scattered symbols. Our results calculated from the multi-source idealgas model 
are shown with the curves. 
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FIG. 2: P T dependence of vj for charge hadrons in Pb-Pb collisions at -\Js NN = 2.76 TeV. Experimental data taken from the 
ALICE Collaboration [l]] are shown with the scattered symbols. Our results calculated from the multi-source model idealgas 
are shown with the curves. 
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FIG. 3: Centrality and the impact parameter b dependence of the square of the expansion factor k 2 . The symbols represent the 
parameter values used in Figs. 1-2. The lines are a fitted results. 



